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Abstract 

We present a global fit to HERA data on the reduced cross section measured in electron- 
proton collisions in the region of small Bjorken-x: x < xq = 10~ 2 and moderate to high values 
of the virtuality Q 2 < Qm ax = 650 GeV 2 . The main dynamical ingredients in the fits are two 
recently proposed improved BK equations for the description of the small-x evolution of the 
dipole scattering amplitude [1,2]. These two new equations provide an all-order resummation 
of double collinear logarithms that arise beyond leading logarithmic accuracy. We show that 
a very good description of data is possible in both cases, provided the parent dipole or 
smallest dipole prescriptions are employed for the running of the coupling. 


1 Introduction 

The Balitsky-Kovchegov (BK) equation [3,4] for the small-x evolution of color dipoles provides 
the main practical tool for phenomenological studies of non-linear, saturation effects in available 
experimental data. The BK equation corresponds to the large-fV ( , limit of the B-JIMWLK [4-7] 
hierarchy of coupled non-linear evolution equations. These renormalisation group equations are 
one main pillar of the Color Glass Condensate (CGC) effective theory for the study of QCD 
scattering process in the high-energy regime. 

At leading logarithmic accuracy (LL) the BK equation resums large logarithms ~ a s Y that 
arise from the successive emission of small -x gluons. At this degree of accuracy the rapidity 
variable can be defined as Y = ln(xo/x), where x is the Bjorken variable and Xq a reference, initial 
scale. Additionally, the BK equation also accounts for non-linear, recombination effects that tame 
or saturate the growth of gluon densities with decreasing x. Over the last years great efforts have 
been devoted to systematically improve the dynamical input of the BK and B-JIMWLK equations 
beyond leading logarithmic accuracy. One big step in this direction was the calculation of running 
coupling corrections to the BK and B-JIMWLK evolution kernels performed in [8,9]. There, 
the evolution kernel to running coupling accuracy was obtained via an all-order resummation of 
a subset of next-to-leading corrections (NLO). More recently, full NLO corrections to the BK 
equation have been calculated in [10]. Progress in the refinement of the theoretical formulation of 
the CGC is not limited to the study of evolution equations. Thus, full or partial next-to-leading 
corrections to factorisation theorems for particle production in deep inelastic scattering (DIS) and 
hadronic collisions have been recently presented in [11-16]. 
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Despite the notable progress briefly reported above, it turns out that the perturbative ex¬ 
pansion devised in these works shows an unstable behaviour. Phenomenological studies on the 
effect of partial [17,18] and full [19] NLO corrections in single inclusive particle production in 
proton-proton and proton-nucleus collisions found that those corrections become very large in cer¬ 
tain regions of phase space, large transverse momentum and small rapidities, overcoming the LO 
contribution and even yielding negative values of the cross sections studied in those works. These 
findings on the instability of NLO corrections have been confirmed recently in a first numerical 
study of the solutions of the full BK equation at NLO presented in [20]. There, it was found 
that the solutions of the BK equation at NLO accuracy turn negative for a large class of initial 
conditions relevant for phenomenological applications. 

The lack of convergence and unstable behaviour of the perturbative expansion of high-energy 
evolution has been identified as due to the presence of large transverse logarithms that arise 
beyond the leading logarithmic approximation [20]. A related problem -large negative contri¬ 
butions arising from transverse logarithms- was identified after the NLO calculation of BFKL 
evolution [21-23], the linear limit of the BK equation. 

Although their precise derivation is rather technical, the problematic nature of these logarithms 
can be understood on general physical grounds: The derivation of the BFKL and BK equation 
rely on the strong ordering of successively emitted gluons. This separation of time scales paves the 
way for the factorisation of high-energy evolution: new emissions do not alter the kinematics of 
previously emitted gluons and can be regarded as independent, i.e their emission factorises. Given 
that the lifetime of an emitted gluon is r = k + /2k“j_, this ordering is achieved by imposing strong 
ordering of the plus light-cone component of successive emissions q + kf ... fc+, provided 
that their transverse momenta are of the same order 1/r ^ Qo ~ k±\ ~ ... k± n . Here 1/r and Q o 
are the characteristic transverse momentum scales of the projectile and target, respectively. In 
the dipole model of DIS they can be related to the transverse dipole size, r ~ l/Q, where Q is the 
photon virtuality, and the saturation scale of the target Qq. However, the time order condition 
at the basis of high-energy factorisation is violated by emissions of large transverse momentum, 
which contribution is parametrically controlled by the the logarithmic extent of the transverse 
phase space p = ln(l/rQo); and become as large as the leading logarithms when Y ~ p. Formally, 
these contributions arise at NLO accuracy as double-logarithmic corrections ~ a s Yp and a s p 2 , 
globally referred to as double collincar logs. Clearly, their importance increases for large Q-values, 
i.e in the collinear limit, and their adequate treatment is crucial to match properly with DGLAP 
evolution. 

An alternative approach to standard order-by-order calculations consists in rearranging the 
perturbative expansion via a resummation of large double logarithmic corrections to all orders. 
Such has been the strategy followed in two recent works: [1] and [2], The net effect of these 
resummations is to limit the longitudinal and transverse phase space in order to ensure the time 
ordering condition that lies at the basis of the BK and BFKL equations, thereby restoring the 
stability of the perturbative series. As we shall see, a direct consequence of this reduced phase 
space is the further slow down of the evolution speed. The work of [1] relies in the use of the 
Mellin space approach and results in the modified BK equation Eq. (2.7). This equation explicitly 
includes a kinematic constraint in the evolution kernel in the form of theta function that reduces 
the transverse phase space and is non-local in rapidity. We shall refer to it as kinematically 
corrected BK equation, KC-rcBK. In turn, the collinearly improved BK equation presented in [2], 
Eq. (2.9), relies in a diagrammatical calculation, was fully derived in coordinate space and is local 
in rapidity. The locality of the improved BK equation is achieved via an analytic continuation of 
the physical scattering amplitude outside the physics range Y > p. We shall refer to this equation 
as DLA-rcBK in what follows. 

Here we aim at exploring the compatibility of those two collinearly improved BK equations 
with data on the reduced cross section in e+p collisions at small values of Bjorxen-x measured 
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at HERA by the HI and ZEUS collaborations [24], The role on non-linear corrections in small-x 
QCD evolution in the interpretation of HERA dada has been thoroughly investigated in previous 
works [25-29] . The main conclusion extracted from these works is that the BK equation including 
only running coupling corrections to the evolution kernel provides a very good description of 
all available experimental data at small Bjorken-x. Further, the precise quantitative information 
extracted from these fits in the form of parametrizations of the dipole-proton scattering amplitude 
has become an essential tool to calibrate the physics expectations and to analyse data from 
other experimental programs, most notably from the proton-proton, proton-nucleus and nucleus 
collisions performed at RHIC and the LHC (see e.g. [18,30]). 

During the completion of this work a similar analysis of HERA data based on the DLA-rcBK 
equation Eq. (2.9) has been presented in [31]. The fits to HERA data presented in that work rely 
on a very similar set up to the one followed here, the main methodological differences being the 
choice of initial conditions and the prescription employed for the running of the strong coupling. 
Further, the work presented in [31] also explores the role of single transverse logarithms ~ a s p in 
the description of data. Nonetheless, the results of this work in the form of parametrisations of 
the dipole-proton scattering amplitude have already been employed in a previous publication [32] 
for the calculation of the neutrino-nucleon cross section at ultra-high energies. 

This work is structured as follows: In section 2 we discuss the basic elements of the collincarly 
improved BK equations proposed in [1] and [2]. We shall refer to these two equations as kinemati¬ 
cally corrected running coupling BK equation (KC-rcBK) and DLA running coupling BK equation 
(DLA-rcBK) respectively. In section 3 we review briefly the dipole model of deep-inelastic scatter¬ 
ing and describe the numerical set up for the fits, including the description of the free parameters 
to be fitted to experimental data. The only novelty with respect to the procedure employed in 
previous works [25-27] is the use of a new family of initial conditions for the BK evolution that 
we shall refer to as pre-scaling initial conditions. Finally in sections 4 and 5 we present the results 
of the fits and the conclusions. 


2 Collinearly improved BK equations 


In this section we briefly review the basic elements of the BK equation including both running 
coupling and double logarithmic corrections. Before delving into details, it is important to recall 
that running coupling corrections and double logarithmic corrections arise from the resummation 
of two different subsets of next-to-leading terms. Therefore, they can be treated as independent 
from each other and straightforwardly combined in a single equation by, for instance, replacing 
the LL kernel in the improved collincar equations by the running coupling kernel. This is exactly 
the procedure to be followed in this work. 

The scattering matrix of a colourless quark-antiquark dipole propagating through the gluon 
held of a target hadron reads: 

S(x 0 ,Xi;K) = ^r(tr {U(x 0 )f/ t (xi)}) y = 5 0 i ; r (2.1) 

where x 0j (i) are the transverse coordinates of the quark and antiquark respectively. Under the 
eikonal approximation the propagation of each individiual right-moving parton through the hadron 
target is accounted for time ordered Wilson lines: 


f/(x) = P exp 




( 2 . 2 ) 


where A + denotes the gluon held of the target hadron. The average in Eq. (2.1) is performed 
over the target gluon held configurations at a given rapidity Y = ln(xo/x). We have used the last 
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line of Eq. (2.1) to introduce the reduced notation for the average of the dipole scattering matrix 
that we shall employ hereafter. The BK equation provides the rapidity evolution of the dipole 
scattering matrix. It reads 

■M-012 [«So2 ; y <Si2,Y — <S 0 i ; y] j (2-3) 


dSpi-Y _ r d 2 x 2 

dY ~ J 2vr 


where .A/1 012 is the evolution kernel. The BK equation at leading logarithmic accuracy resums large 
logarithmic corrections of the form ~ ( a s Y) n to all orders. It accounts for the multiple emission of 
soft gluons and also for the possibility of gluon recombination via the non-linear term in the right 
hand-side of Eq. (2.3). The BK equation assumes that the average of the product of two dipole 
scattering matrices over the held configurations of the target factorizes: (S xz S zy )Y — > <S xz; y<S' zy; y, 
i.e. it is derived in the mean held limit. We shall rely on the translational invariant approximation, 
equivalent to the assumption that de dipole amplitude only depends on its relative transverse size. 
We tus introduces the variables r 0 = |x 0 — xjJ, rq = |x x — x 2 | and r 2 = |x 0 — x 2 |. 

The calculation of running coupling corrections to the original LL kernel was performed in [8,9]. 
There, two different prescriptions were proposed for the kernel of the BK equation at running 
coupling accuracy. It was shown in [33] that Balitsky’s prescription [9] minimizes the role of 
higher-order, conformal corrections, suggesting that it may be better suited for phenomenologi¬ 
cal applications. In particular, Balitsky’s prescription Eq. (2.4) yields a slower evolution speed 
than other possible schemes explored in the literature, like the Kovchegov-Weigert [8] one or the 
smallest dipole size prescription, where the scale for the running coupling is set by the smallest 
of the transverse dipole sizes involved in one evolution step: r 0 , parent dipole and r 12 , daughter 
dipoles. This feature was crucial for the very good description of previous AAMQS fits to HERA 
data. However, it remains to be clarified which running coupling scheme is theoretically better 
motivated once other dynamical effects, like the double logarithmic corrections discussed here, are 
also incorporated to the BK equation. Under Balitsky’s prescription the running coupling kernel 
reads 


M 


Bal 

012 


a s (rl)N c 

7T 


, I ( <r\) \ 1 ( as{rp _ Y 

r\rl + rl\a s (rl) ) + r 2 W?) )_ 


Bal. (2.4) 


Here we will consider two other possibilities for the running coupling kernel, namely the parent 
dipole prescription: 
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(2.5) 


and the smallest dipole prescription: 
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012 
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7T 
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rrr, 


1 ' 2 


with r min = min{r 0 , rq, r 2 } 


SD. 


( 2 . 6 ) 


Although the ansatzs Eq. (2.5) and Eq. (2.6) do not follow from any strict diagrammatical calcu¬ 
lation, it was shown in [33] that the parent dipole description leads to solutions very similar to 
those obtained under the prescription derived by Kovchegov and Weigert [8]. In particular, the 
parent dipole prescription leads to significantly faster evolution speed that Balitsky’s prescription. 
We shall use it here, rather than the full Kovchegov-Weigert kernel, due to its relative simplicity in 
the numerical evaluation. In turn, the smallest dipole prescription is motivated by the expectation 
that the scale for the running coupling should be given by the hardest momentum scale in the 
process. 
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2.1 Kinematically corrected BK equation (KC-rcBK) 

A kinematically improved version of the BK equation consistent at high, but finite, energies has 
been proposed in [1], It reads: 

QY 5 = j ^ A4 0 12 — A012) [5 0 2 ; r-Aoi2 ‘5l2;Y-A 0 i2 — £oi ; y] • (2.7) 

The reduction of phase space for the evolution is made explicit by the presence of the theta 
function in Eq. (2.3), which constrains gluon emission to some bounded domain in the x 2 -plane. 
This is in clear contrast with the LL BK equation or the BK equation with only running coupling 
corrections, where emission of gluons in all the transverse plane are allowed. The parameter that 
controls the extent of kinematic corrections is given by: 

A 012 = max jo, In j witl1 ^012 = min (ri, r 2 } . (2.8) 

There is actually some freedom in the definition of Z 0 i 2 introduced in Eq. (2.8). It should satisfy 
Z 012 ~ rq ~ r 2 in the regime r 0 <C rq ~ r 2 . This freedom in the choice of A 0 i 2 should be considered 
as a resummation scheme ambiguity associated with the kinematic constraint. In what follows 
we shall adopt the definition presented in the r.h.s of Eq. (2.8); we have checked that changes in 
this prescription do not alter significantly the results of the fits. Other important feature of the 
kinematically corrected BK equation Eq. (2.8) is that the scattering amplitude of the two newly 
created dipoles after one evolution step are evaluated at a delayed rapidity Y — A 012 . This rapidity 
veto has previously discussed in the literature as a main part of the NLO or energy-momentum 
corrections to the BK equation [28,34], As we shall see in Section 4 both effects tend to decrease 
the evolution speed, i.e. to slow down the growth of the saturation scale with decreasing Bjorken- 
x, but they also modify more exclusive features of the unintegrated gluon distribution. Another 
subtle point that arises in the definition of the high-energy factorisation scheme once the kinematic 
corrections are taken into account is related to the very definition of the evolution variable, i.e. 
the rapidity variable. We will not delve into the details here (see the extended discussion in [1]). 
It is however important to recall that the usual choice Y = In (k~jt/k^) = \n(xo/x) can be modified 
by finite corrections, leading to different evolution equations at NLO accuracy and beyond. This 
freedom to choose the evolution variable is related to freedom to choose a reference energy scale 
in the BFKL formalism. In practice we will deal with this ambiguity through the use of pre- 
asymptotic initial conditions, see 3. The introduction of an arbitatry rapidity shift AK 0 as another 
free fit parameter copes effectively with a possible redefinition of the rapidity variable. 


2.2 BK equation at Double Logarithmic Accuracy (DLA-rcBK) 


A BK equation that resums double collinear logs to all orders has been recently proposed in [2], 
It reads: 


<9<Soi ; y 
d Y 


d 2 x 2 
2vr 


DLA 

012/V 012 


At 012 A, 


*Si 2; y 


5, 


01;Y 


(2.9) 


It should be noted that the object evolved in Eq. (2.9), is not the physical dipole scattering matrix 
5oi ; y, but rather a related function 5 0 i ; y. Both functions are expected to coincide in the physical 
range Y > p. For general positive values of Y and p the relation between these two quantities is 
given by 

A{Y,p)= f P dpJfrp-pJAfrpi) ( 2 . 10 ) 

Jo 
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with 


f(Y = 0 ,p) = S(p) - y^J 1 (2 V / ^?), (2.11) 

where Ji is the Bessel function of the first kind. The function A(Y, p) in Eq. (2.10) is defined as 
(1 — iS xy; y) = r 2 QlA xy -Y, with p = \i\{l/r 2 Ql) and Q o some initial scale (its precise definition 
is given in section 3). The introduction of the auxiliary function <S 01 .y obeys the purpose of 
obtaining a collinearly improved equation that, contrary to Eq. (2.7), is local in rapidity. This 
implies performing an analytic continuation of the dipole scattering matrix outside the physical 
range Y > p. Eq. (2.9) exhibits two important features. First, it is local in rapidity. Second, 
the resummation of double logarithmic corrections results in just a modification of the evolution 
kernel by the factor 


Kgff = with p'= 7 ln M/ r »> In (r?/rg). (2.12) 

It is important to note that the resummation affects both the evolution kernel and the initial con¬ 
ditions for the evolution. The latter, A, are obtained applying Eq. (2.10) to the initial conditions 
defined in Eqs. (3.5) - (3.6) in Section 3. Finally, the derivation of [2] treats the strong coupling as 
a fixed parameter. In order to be consistent with the degree of accuracy of the evolution kernel(s) 
Eq. (2.4) and Eq. (2.5) we shall let the coupling run in Eqns. (2.10)-(2.12) at the scale p = 2/r 
or p = 2/min(r, rq, r 2 ). As we shall discuss in section 4, it turns out that the results of the fits 
are very little sensitive to this choice. 


3 Set up 


3.1 Dipole model of DIS 

In this section we briefly review the main ingredients needed for the calculation of the inclusive 
and longitudinal DIS structure functions, which was extensively discussed previous papers, see 
e.g. [25,26]. Neglecting the contribution from Z boson exchange, the reduced cross section can be 
expressed in terms of the inclusive, F 2 , and longitudinal, F L , structure functions: 

a r (y,x,Q 2 ) = F 2 (x, Q 2 ) - - — J _ F L (x, Q 2 ), (3.1) 


where y = Q 2 /(s x) is the inelasticity variable and a/s the center of mass collision energy. In turn, 
at x <C 1, the inclusive and longitudinal structure functions can be expressed as 


F2(x,Q 2 ) 


Q 2 

4 7T CY em 


(a T + cr L ) , 


Fl(x,Q 2 ) 


Q 2 

TY -Tl ■ 

4 7T Ot e m 


(3.2) 


Here ^ stands for the virtual photon-proton cross section for transverse (T) and longitudinal 
(L) polarization of the virtual photon. In the dipole model, valid at high energies or small x, one 
writes [35,36]: 

<x T ,l( x ’Q 2 ) = ^2 J o dz J d 2 r \^T,L( e f’ m fA,Q 2 ,r)\ 2 <x qq (r,x),, (3.3) 


where L ^ 1C light-cone wave function for a virtual photon to fluctuate into a quark-antiquark 
dipole of quark flavor /. Note that ^ TL only depends on the quark flavor / through the quark 
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mass rrif, and electric charge e.j (see e.g. [37] for explicit expressions to lowest order in a em ). 
According to the optical theorem, the dipole cross section a qq ( r, x) is given by the integral over 
impact parameter of the (imaginary part of) dipole-hadron scattering amplitude. Further, under 
the approximation that the dipole scattering amplitude is independent on the impact parameter 
of the collision one gets 


cr qq (r,x) = 2 J d 2 bj\f{x,r,b ) = a 0 JV(x,r) (3.4) 

where 07 has the meaning of (half) the average transverse area of the quark distribution in the 
transverse plane and will be one of the free parameters in the fit. Following [26] we have consider 
a variable number of active flavours up to Nf— 5, with a current quark mass m u ,d, s ,c,b = 0.05, 0.05, 
0.14, 1.3 and 4.5 GeV respectively. It turns out that the fits results are very little sensitive to 
modifications of the particular choice made for the current mass of light flavours. 


3.2 Pre-asymptotic initial conditions 

In order to solve the BK equations under the different evolution schemes discussed in section 2 
one has first to specify the initial conditions for the evolution at the highest value of Bjorken- 
x included in the data fitting set, Xq = 0.01. The free parameters in the AAMQS fits mainly 
correspond to the free parameters in the initial conditions for the evolution, as follows: 


A f(r, Y = 0) 


1 — exp 



In 



(3.5) 


Q o and 7 are two free parameters to be fitted to experimental data. 7 is a dimensionless parameter 
that controls the steepness of the unintegrated gluon distribution at moderate and high transverse 
momenta. I 11 the original AAMSQ fits, Q 0 plays the role of the initial saturation scale at the 
highest value of Bjorken-x included in the fitting data set, x 0 = 1CT 2 or, equivalently at Y — 0 
with Y = ln(x 0 /x). However, it turns out that it is not possible to obtain good fits for the 
collinearly improved equations discussed here using the two-parameter family of initial conditions 
provided by Eq. (3.5). I 11 order to allow a larger freedom for the functional forms of the initial 
conditions we shall introduce a new family of pre-asymptotic initial conditions as follows: We first 
solve the BK equation using Eq. (3.5) as the initial conditions for the evolution and take the dipole 
scattering amplitude at x = x 0 as the solution of the BK equation at rapidity AK 0 : 


A f(r, x 0 = 0 . 01 ) = A f(r, A Y 0 ). 


(3.6) 


That is, we allow the evolution to run for some rapidity interval AYo before comparing to ex¬ 
perimental data. In other words, we use the evolution itself to generate the initial conditions for 
further evolution at Y > AF 0 or, equivalently for x < x 0 . Thus, the dipole scattering amplitude 
in the physical region x < Xq is given by 


A7(r, x < xq) = A f(r, AT'o + ln(x 0 /x)), 


(3.7) 


where AK 0 quantifies the amount of pre-evolution allowed before comparing to experimental data 
and is another fitting parameter. Several comments are in order: 

First, this set up should just be regarded as a mathematical procedure to generate a family 
of initial conditions which are a solution of the BK equation itself. Therefore, for values AYo > 0 
neither the fitting parameters Q 0 and 7 nor the solutions of the BK equation for Y < AYo afford a 
clear, straightforward physical interpretation. Thus the scale Qo in Eq. (3.5) is just an auxiliary 
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scale to generate the physical initial conditions at x 0 = 1CT 2 . The physically meaningful object 
is the dipole scattering amplitude for x < x 0 or Y > AV 0 as defined in Eq. (3.6) and Eq. (3.7). 
Analogously, the physically meaningful saturation scale is defined via the condition 

Af(r — 1/Q s (x), x) — 0.5 (3.8) 


for x < xq. 

Second, similar initial conditions have been previously used in the literature in the limiting case 
of AYo —>> oo. These are referred to as scaling initial conditions in [25,28]. It is well known that 
the non-linear character of the BK equation leads to scaling solutions, i.e independent of the initial 
conditions at asymptotically large rapidities. In that limit, the solutions depend only on a single 
variable, the scaling variable r = rQ s (x), hence effectively becoming a become a one-parameter 
family of possible initial conditions. This property of the BK equation has been studied in depth 
in relation to the observed phenomenon of geometric scaling featured by experimental data on 
e+p collisions from HERA. Here we take an intermediate approach where we allow for a finite 
amount of evolution AY 0 before comparing to data. In this intermediate limit the initial conditions 
generated at x = x 0 according to Eq. (3.6) are still sensitive to the initial fitting parameters Q 0 
and 7 . 

Finally, the pre-evolution interval AYo helps the convergence of the solution of the DLA equa¬ 
tion to its regime of physical applicability, defined by the condition Y > p = ln(l/r 2 QQ) and, in a 
related way, this approach copes with the uncertainty in the definition of the evolution variable, 
i.e the rapidity Y, in the evolution scheme including kinematic corrections. 


3.3 Regularisation of the infrared dynamics 

The solution of the BK equation in either of the evolution schemes described above implies the 
evaluation of the strong coupling at arbitrarily large values of the dipole size (small gluon momen¬ 
tum), and a regularization prescription to avoid the Landau pole becomes necessary. Following 
previous AAMQS works, for small dipole sizes r < r/ r we shall evaluate the running coupling 
according to the usual one-loop QCD expression: 



12vr 

(nw.-awin^j 


(3.9) 


The number of active flavors Nf in Eq. (3.9) should is set to the number of quark flavors lighter than 
the momentum scale associated with the scale r 2 at which the coupling is evaluated p 2 = 4C 2 /r 2 . 
The setup of this variable flavor scheme is completed by matching the branches of the coupling 
with adjacent Nf at the scale corresponding to the quark masses r 2 = 4C 2 /m 2 . For the 1-loop 
accuracy at which the coupling Eq. (3.9) is evaluated, the matching condition is simply given by 


oc s ,N f -i(r *) = ot SjNf (n) . (3.10) 

With only three active flavours one gets A q C d = 0.241 GeV, such that a s (M z ) = 0.1176, with 
M z the mass of the Z boson. In turn, for larger sizes, r > ry r , we freeze the coupling to the fixed 
value a s (rf r ) = af r . We shall use two different values aj r = 0.7 and 1. The fudge factor C under 
the logarithm in Eq. (3.9) will be one of the free parameters in the fit. It reflects the uncertainty in 
the Fourier transform from momentum space, where the original calculation of a s Nf corrections 
was performed, to coordinate space. We shall extend the fits up to values of the photon virtuality 
Q 2 ~ 650 GeV 2 . 
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Qmax (GeV^) 

Evolution 

scheme 

a fr 

^(GeV 2 ) 

AK 0 

(To (mb) 

7 

C 

X 2 /d.o.f. 

50 

Bal 

0.7 

0.210 

0 

24.03 

1.141 
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0.881 

Bal 

0.7 

0.182 

0 

24.06 

1.071 

1.36 

1.119 

Bal 

1 

0.213 

0 

24.02 

1.147 

1.147 

0.875 

Bal 

1 

0.197 

0 

25.00 

1.101 

1.39 

0.924 

650 

Bal 

0.7 

0.216 

0 

24.16 

1.160 

1.461 

0.928 

Bal 

0.7 

0.209 

0 

23.14 

1.124 

1.250 

1.08 

Bal 

1 

0.218 

0 

23.94 

1.158 

1.232 

0.869 

Bal 

1 

0.207 

0 

23.14 

1.122 

1.12 

1.119 


Table 1: Fit parameters for rcBK evolution with the running coupling kernel evaluated according 
to Balitsky’s prescription Eq. (2.4). 
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PD 

1 
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1.13 

SD 

0.7 

5.38-10" 2 
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1.166 
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Table 2: Fit parameters for DLA-rcBK evolution with the running coupling kernel evaluated 
according to the Bal, PD and SD schemes. 
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scheme 

°fr 

Q*(GeVg 

AKo 

<T 0 (mb) 

7 

C 

X 2 / d.o.f. 

50 

PD 

0.7 

4.65-10" 2 

4.98 

22.33 

1.001 

3.662 

1.081 

Bal 

0.7 

1.66-10" 2 

5.53 

23.75 

0.845 

0.869 

1.332 

PD 

1 

3.35-10" 2 

6.60 

22.23 

1.032 

3.806 

1.108 

Bal 

1 

2.07-10" 2 

6.15 

22.85 

0.8542 

0.675 

1.246 

650 

PD 

0.7 

5.22-10" 2 

5.03 

22.79 

1.058 

3.67 

1.372 

PD 

1 

0.33T0~ 2 

6.14 

22.19 

1.040 

3.74 

1.280 

SD 

0.7 

0.101 

3.24 

21.55 

1.101 

2.52 

1.21 


Table 3: Fit parameters for KC-rcBK evolution with the running coupling kernel evaluated ac¬ 
cording to the Bal, PD and SD schemes. 
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4 Results 


With the set up described in the previous section it is possible to obtain good fits to data on 
the reduced cross sections with either of the evolution schemes rcBK, KC-BK and DLA-BK. The 
corresponding fit parameters and the y 2 jd.o.f are shown in Table 1-3, whereas a comparison 
between the fit results and experimental data for some selected Q 2 bins is shown in Fig. 2. 

Before addressing specific aspects of the fit results, it should be noted that all successful fits 
(those with y 2 /d.o.f. < 1.25) yield a very similar value of the dipole scattering amplitude in the 
region of Bjorken-x and Q 2 constrained by experimental data. This is not an unexpected result 
since, once the calculation framework is fixed -the dipole model in our case- the simultaneous 
description of the x and Q 2 -dependence of data leaves little room for variations of the optimal 
dipole amplitude. Indeed, the only fit parameter that does not relate directly to the behaviour of 
the dipole amplitude, either of its initial condition or its evolution, is the normalisation parameter 
(Jo- We observe that the variations of the dipole amplitude among different fits are compensated by 
changes of the same order in the corresponding values of the 07 parameter, ~ 10%. The mentioned 
compatibility among different evolution schemes is clearly seen in Fig. 1, left plot, which shows the 
dipole scattering amplitude resulting from four different fits at the highest x -value included in the 
fitted data set as a function of the dipole transverse size. The four dipole amplitudes represented 
in that plot deviate little from each other in the region of large and moderate dipole sizes. Using 
the ballpark estimate that relates the photon virtuality and the dipole size as Q ~ 2/r, fits up to 
Qmax = 500 GeV 2 (there is a single point in the fitting set with Q 2 > 500 GeV 2 ) constrain the 
region of dipole sizes r > 1CT 1 GeV -1 . It is precisely outside this region, i.e for r < 10 -1 GeV -1 , 
that differences between the different parametrisation of the dipole amplitude start growing above 
10 % and higher. 

Fig. 1, right plot, shows the value of the saturation scale, obtained according to the definition 
given in Eq. (3.8), as a function of the rapidity variable Y — \w(;xq/x) for four different fits. Again, 
and in agreement with the previous discussion on the similarity of the dipole amplitudes, we 
observe that different fitting schemes lead to a very similar value for the saturation scale. In 
particular they all yield an initial value Q 2 (x 0 = 0.01) ~ 0.2 GeV 2 . This value is in very good 
agreement with the one obtained in previous fit to HERA data. Differences in the rate of growth of 
Q 2 (x ) with decreasing x are larger than those of the the dipole amplitude itself, with the rcBK-Bal 
scheme yielding the slowest evolution. 

For reference we have first performed fits using the BK equations with only running coupling 
corrections. The fit parameters are shown in Table 1. In this case it was possible to obtain good 
fits using only Balitsky’s scheme for the running of the coupling. The other two schemes, PD and 
SD, lead to a too fast evolution speed that can not be reconciled with data. The fits are stable 
when changing the maximum virtuality in the fitted data set from 50 to 650 GeV 2 . The best fits 
are obtained when AVo is fixed to 0 and no restriction is imposed on the parameter 7 , which tends 
to acquire values ~ 1.14 -7 1.15. It turns out that these range of 7 values lead to a non positive 
definite fourier transform of the dipole amplitude and, hence, to an unphysical unintegrated gluon 
distribution. In order to guarantee the physicality of the related gluon distribution we impose 
the constraint 7 < 1.125. We checked that below this limiting value the fourier transform of the 
initial conditions Eq. (3.5) is positive definite. This feature is preserved by the evolution, as it 
tends to decrease the effective anomalous dimension of the solutions down to its asymptotic value 
7 ~ 0.85 [38,39]. This additional constraint leads to slightly higher y 2 /d.o.f, but still or order 
0(1) and stable against increasing Qm ax - It should be noted that the fits presented here vary 
slightly with respect to analogous fits presented in [25,26]. The main reason for this change of 
the fits results is the improvement of the numerical accuracy of the code and, mostly, to a more 
refined search strategy in the fitting algorithm that allows to find a larger amount of local minima. 

Similarly, a very good description of data can be obtained using the DLA-rcBK equation 
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Eq. (2.9) for the evolution of the dipole amplitude. The fit parameters are shown in Table 2. In 
this case, Balitsky’s scheme for the running coupling yields the worst fits to data, y 2 /d.o.f > 1.5. 
In turn, it is possible to get very good, stable with increasing Q 2 max fits with he other two running 
coupling schemes, PD and SD. This can be understood as due to the lower evolution speed yielded 
obtained within Balitsky’s prescription with respect to the PD and SD schemes. The evolution 
speed is further reduced when the additional phase space restrictions entailed by the collinear 
resummation of Eq. (2.9) are also incorporated in the evolution kernel. The most remarkable 
feature of these fits is that they require large pre-evolution periods: 2.9 < AYq < 9.1. This 
brings the physical solutions deeper into the scaling regime, thus speeding up the evolution and 
compensating the speed reduction induced by DLA effects. The fudge factor C takes values larger 
than for rcBK fits: C ~ 2.5 -j- 3.8. We have checked that, as suggested in [31], the effects of the 
DLA resummation on the initial conditions can be reabsorbed in to a slight reshuffling of the fit 
parameters. Nonetheless, all the fits described here have been obtained including the resummation 
effects on the initial condition, i.e. applying Eq. (2.10) to the initial conditions defined in Eqs. 
(3.5) - (3.6). The particular choice for the scale setting the running of the coupling in Eqns. 
2.10-2.12 did not have any significant impact on the fit results. It should also be noted that the 
very small values of the initial scale Qo, only slightly larger than A qcd, do not carry any special 
physical meaning. Rather, the physically meaning full dipole amplitudes and the corresponding 
saturation scales are very similar to those obtained with rcBK evolution only, as already discussed 
at the beginning of this section and illustrated in Fig. 1. 

Finally, fits using the kinematically corrected BK equation KC-rcBK Eq. (2.7) also provide a 
very good description of the data. The fit parameters are shown in Table 3. The discussion of the 
results here runs close to the one for the DLA-rcBK fits. Again, the use Balitsky’s scheme yields 
reasonable fits only for Q^ nax = 50 GeV 2 , turning unstable for larger photon virtualities. In turn 
the use of PD or SD schemes for the running coupling kernel results in good, stable fits. However, 
and contrary to the rcBK and DLA-rcBK approaches, the quality of the fits worsens slightly 
when Qlnax is increased, revealing some tension when the fits are pushed towards the collinear 
limit. Again, KC-rcBK fits require a considerable amount of pre-evolution before comparing to 
experimental data: < 3.2 < AK 0 < 6.6, while the values obtained for the fudge factor are very 
similar to those for DLA-rcBK fits: C ~ 2.5 3.7. 



Left: Dipole scattering amplitude versus dipole size at x — 0.01 for four of the fit parametri- 
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sations. Right: Saturation scale corresponding to the dipole amplitudes represented in the right 
plot as a function of Bjorken-x 



Figure 2: Comparison between data on the reduced cross section (full circles) and the fit results 
for three different evolution schemes: rcBK-Bal (open circles), DLA-rcBK-PD (open squares) and 
KC-rcBK-SD (open triangles). The three fit results correspond to «/ r = 0.7. The lower panels in 
each plot show the ratio theory/data together with the experimental uncertainty. 


5 Conclusions 

To conclude, in this work we have shown that the two collinearly improved BK equations presented 
in [1, 2] provide a good description of HERA data on the reduced cross section measured in 
electron-proton collisions in the small Bjorken-x region, x < 0.01, and for moderate to large 
virtualities, Q 2 < 650 GeV 2 . This has been possible through the introduction of a new set of 
pre-asymptotic initial conditions. This exercise alone does not suffice to distinguish which of 
the two approaches studied here is better suited to describe the high-energy evolution of color 
dipoles. Rather, additional theoretical effort is needed in order to determine the right running 
coupling evolution scheme and also to correct the factorisation theorems for particle production 
to the same degree of accuracy. Also, the phenomenological analysis of more exclusive processes is 
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needed in order to shed more light on the compatibility of the different resummation procedures 
with experimental data. 
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